Constraints on Axial Two-Body Currents from Solar Neutrino Data 
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We briefly review recent calculations of neutrino deuteron cross sections within the effective field 
theory and traditional potential model approaches. We summarize recent efforts to determine the 
counter term describing axial two-body currents, Lia, in the effective field theory approach. We 
determine the counter term directly from the solar neutrino data and find several, slightly different, 
ranges of Lia under different sets of assumptions. Our most conservative fit value with the largest 
uncertainty is Lia ~ 4.5^2 f m3 - We show that the contribution of the uncertainty of Lia to the 
analysis and interpretation of the solar neutrino data measured at the Sudbury Neutrino Observatory 
is significantly less than the uncertainty coming from the lack of having a better knowledge of #13 . 
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A significant amount of theoretical work was recently 
directed towards the calculation of neutrino capture on 
deuteron. Some of these efforts to describe this process 
utilized the effective field theory approach. In this ap- 
proach nonlocal interactions at short distances are rep- 
resented by effective local interactions in a derivative ex- 
pansion. Since the effect of a given operator on low- 
energy physics is inversely proportional to its dimension, 
an effective theory valid at low energies can be written 
down by retaining operators up to a given dimension. 
The coefficients of these operators are then needed to be 
fixed either directly by the data or can be fitted to the 
results of calculations carried out using more traditional 
approaches. 

For nucleon-nucleon interactions it was shown that one 
can introduce a well-defined power counting [|. In this 
method one needs to introduce a single coefficient, com- 
monly called L\Ai to parameterize the unknown isovector 
axial two-body current which dominates the uncertain- 
ties of all neutrino-deuteron interactions. Using an ef- 
fective theory without pions 0] such a calculation was 
carried out in Ref . .3] . These authors found that the ra- 
tio of charged- to neutral-current was fairly insensitive to 
this counter term. To test the convergence of the results 
in Ref. [3 Butler, Chen, and Kong also calculated next- 
order corrections and found that no new parameters need 
to be introduced Q|. An alternative formulation of the 
effective field theory approach using heavy-baryon chiral 
perturbation theory was given in Ref. [{| . 

The cross section for neutrino absorption on deuterium 
was first calculated in Refs. [|| and utilizing an effec- 
tive range approximation to describe the nuclear interac- 
tion, using the allowed approximation for the weak oper- 
ators, and assuming that the final two-nucleon state has 
a relative angular momentum of zero. First-forbidden 



contributions to the weak operators were included using 
Siegert's theorem in Ref. [g and using convection cur- 
rent form of the vector operators in Ref. 0- A detailed 
assessment of various approximations in these papers was 
given 111 Refs. [ljj and |lj using various nuclear poten- 
tials. This work was recently updated in Ref. |12| . 

Radiative corrections to the charged-current bre akup 
of the deuteron were calculated b y T owner in Ref. [lj| • 
Beacom and Parke pointed out |I4| inconsistencies in 
Towner's treatment of radiative corrections. This incon- 
sistency was resolved in Ref. and cross section calcu- 
lations u sing more recent values of qa were given in Refs. 
[Isj l and 0| ■ More recently it was shown that radiative 
corrections to the charged-current neutrino-nuclear reac- 
tions with either an electron or a positron in the final 
state are described by a universal function |17j . 

The counter term, Lia, describing the effects of the 
leading weak axial two-body current can be determined 
either by comparing various cross sections calculated us- 
ing the effective field theory approach with those calcu- 
lated using standard potential model approach or with 
experimentally or observationally determined cross sec- 
tions. When the renormalization scale is set to the muon 
mass, dimensional analysis gives a rough estimate of this 
quantity 



6fm d 



(1) 



It should be emphasized that this number depends on 
the renormalization scale and cannot be reliably used at 
lower energies. Using the existing reactor antineutrino- 
deuteron breakup data provides a constraint of |18| 



L\a = 3.6±5.5fm 3 



(2) 
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Helioseismic observation of the pressure-mode oscilla- 
tions of the Sun can be used to put constraints on various 
inputs into the Standard Solar Model, in particular the 
pp fusion cross section. This process has been calculated 
to the fifth order in pionless effective field theory [lflj . 
Neutrino-deuteron and antineutrino-deuteron scattering 
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SNO Allowed Regions for Various Values of l_ 1A 
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FIG. 1: The change in the allowed region of the neutrino parameter space using solar neutrino data measured at SNO as the 
value of Lia changes. In the calculations leading to this figure the neutrino mixing angle #13 is taken to be zero (see text). 
The shaded areas are the 90 % confidence level region. 95 % (solid line), 99 % (log-dashed line), and 99.73 % (dotted-line) 
confidence levels are also shown. 



are computed to the third order in the same approach. 
The value of Lia is not the same in different orders. He- 
lioseismology limits Lia — 7.0±5.9fm 3 in the fifth order 
pfij . Using the expressions given in Ref. 0] this gives 

L XA = 4.8 ± 6.7 fm 3 (3) 

in the third order. A state of the art calculation of the 
pp fusion cross section was given in |2l| where the un- 
certainty in the axial two-body current operator was ad- 
justed to reproduce the measured Gamow- Teller matrix 
element of tritium /3 decay. After performing the trans- 
formation from the fifth- to third-order, the calculation 
of Ref. [2l| indicates a value of 

L\a = 4.2 ± 2.4 fm 3 . (4) 

One can also try to determine the counter term directly 
using the solar neutrino data. Using the Sudbury Neu- 
trino Observatory (SNO) and SuperKamiokande (SK) 
charged-current, neutral current, and elastic scatte ring 
rate data Chen, Heeger, and Robertson (CHR) find [23] 

L\a = 4.0 ± 6.3 fm 3 . (5) 

In order to obtain this result CHR wrote the observed 
rate in terms of an averaged effective cross section and 
a suitably defined response function. In this paper we 
explore the phenomenology associated with the variation 
of L 1A . 

In our calculations we used the neutrino cross sections 
given in Refs. 0] and 0. The radiative corrections are 



taken into account following Ref. [lj|. To calculate ob- 
served solar neutrino rates and spectra we used a co- 
variance approach the details of which are described in 
Ref. [23I • In all calculations to obtain the MSW survival 
probabilities we used the neutrino spectra and solar elec- 
tron density profile given by the Standard Solar Model 
of Bahcall and collaborators |24| . 

The dependence of the extracted neutrino parameters 
on the value of L\a is not very strong. We show how the 
parameter space changes with L\a in Figure 1. In this fig- 
ure to find the allowed regions we fit 34 data points from 
the SNO day-night-spectrum j2^ using the procedure of 
Ref. [2i|. The shaded area is the 90 % confidence level 
region. 95 % (solid line), 99 % (log-dashed line), and 
99.73 % (dotted-line) confidence levels are also shown. 
As Lia changes from —15 fm 3 to 25 fm 3 we note that 
the changes in the shape of the confidence level intervals 
are small. In the calculations leading to this figure we 
took the total 8 B flux to be a free parameter using the 
procedure discussed in Ref. FJU . (Note that even though 
we show the entire parameter space in this figure in the 
rest of this paper we concentrated in the large mixing 
angle region which is preferred by the global analysis). 
We conclude that the uncertainty in Lia can not be a 
significant source of error in the analysis of SNO data. 

In general, both the deuteron breakup cross section 
and the total count rate at SNO are nonlinear in L\a- 
Since L\a is small the charged- and neutral current count 
rates can be linearized by making a first order expansion, 
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FIG. 2: Projection of the global A% 2 function on the param- 
eter Lia- In the calculations leading to this figure #13 is taken 
to be zero. 812 , 8m\ 2 , and the parameter /s (the multiplier 
of the 8 B flux) are varied. The solid line represents the case 
where all these three are unconstrained. The long-dashed line 
is when the 8 B flux is fixed as described in the text, but the 
other two are unconstrained. The dotted line is when #12 
and <5mi 2 are also taken to be the best fit values to the SNO 
energy spectra. 



Count Rate ~ A + B Lia, 



(6) 



as was done by CHR. Both the energy dependence and 
the overall magnitude of the 8 B flux is an input into the 
Standard Solar Model. The energy dependence is rather 
accurately determined by the laboratory measurements 
of the 8 B decay. The overall magnitude is determined 
by the measured rate of the 7 Be(p, 7) 8 B reaction (for 
a review see Ref. 26]). To account for the sensitivity 
of the calculations on the value of the 8 B flux we set 
$( 8 B) = /b < I'ssm( 8 B) and calculate the total rate for 
various values of the parameter fg. Clearly the total 
count rate should be proportional to the value of fg. 
Note that the elastic scattering count rate is independent 
of Lia- Thus allowing fg vary freely cannot be fully 
compensated by changing Lia as we discuss below. 



calculated as a function of L\a- In this figure Ax 2 is 
projected only on one parameter (Lia) so that n — a 
bounds on it are given by A% 2 = n 2 . #13 is assumed to 
be zero. The solid line represents the case in which all 
other parameters (#12, Sm 2 2 , and fg) are unconstrained. 
The best fit value is given by L\a — 4.5 fm 3 . In this case 
L\a is constrained between — 7fm 3 and 23fm 3 at la level. 
Such a wide range is not surprising since the dependence 
of the rate on Lia is small and the effects of the parame- 
ters like #12, Sm 2 2 , and fg are much more dominant. In 
order to obtain a better bound on Lia, we fix fg so that 
the total count rate of SNO [2^ is exactly reproduced 
at the value of Lia which corresponds to the minimum 
X 2 of the fit while the parameters #12 and 8m\ 2 are un- 
constrained. The resulting fit is shown by the dashed 
line. In this case Lia is constrained between —2 fm 3 and 
13 fm 3 at lcr level. The dotted line in this figure repre- 
sents the case where we fix all the parameters except Lia- 
We find the best fit values of 6*12 and 5m 2 2 in a global 
fit using 93 data points from solar and reactor neutrino 
experiments; namely the total rate of the chlorine exper- 
iment (Homestake l2 7| ) . the average rate of the gallium 
experiments (SAGE~|2g, GALLEX GNO wh, 44 
data points from the SK zenith- angle-spectrum |3 J . 34 
data points from the SNO dav-night-spectru m J25l and 
13 data points from the KamLAND spectrum [32|. In ad- 
dition we fix fg so that the total count rate of SNO [2{| 
is exactly reproduced. If we were to exclude SNO data 
from the global analysis and took fg — 1 instead of fixing 
as described above this method would be tantamount to 
treating SNO as an experiment to measure only Lia so 
that the uncertainties in the SNO data would only show 
up as the corresponding uncertainty at Lia- From the 
dotted line Lia is constrained between 2fm 3 and 8fm 3 at 
la level. In a sense this latter range represents the " best 
case" limit on Lia that one can obtain from SNO. It is 
worth emphasizing that the \ 2 minimum is almost the 
same in all these cases. In all cases we obtain a best fit 
value of Lia around 4.5 to 5 fm 3 which is a little larger 
than the value obtained by CHR. These authors use elas- 
tic scattering, charged-current, and neutral-current rates 
separately with effective cross sections. Since we fit the 
solar neutrino day-night spectrum directly by folding dif- 
ferential cross sections, detector response functions, 8 B 
spectrum, and the MSW survival probabilities we need a 
slightly larger Lia- 

In Figure |3| we compare our results with results from 
other analyses. Our results are based on the dashed and 
dotted lines of Figure[21 We calculate la errors by fitting 
a Gaussian of the form 



exp 



r average \ 2 
J 1A 



aL ljf 



(7) 



In Figure |3 we present the quantity A% 2 



X 



Xmin 



to each side of the marginal likelihood expression C = 
exp (— Ax 2 /2) and estimate two standard deviations sep- 
arately for each side j3^| ■ We obtain the error band shown 
in the figure by symmetrizing those errors. 



I_ 1A from Various Analyses 

I l l l l l 
Butler-Chen-Vogel: 3.6 ±5.5 (reactor anti-neutrino) 

Brown-Butler-Guenther: 4.8 ±6.7 (helioseismology ) 
Chen-Heeger-Robertson: 4.0 ±6.3 (CC, NC & ES) 



Balantekin-Yuksel : 4.5 , 2 ( f B , 5m 21 and e 12 are unconstrained) 
Balantekin-Yuksel : 4.5 ±8 (f B is fixed, 8m 21 and 9 12 are unconstrained) 
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FIG. 3: Values of L\a obtained from different analyses. The values labeled this work are calculated using the dashed and the 
dotted lines of Figure 2 as described in the text. Helioseismology limit is from Ref. |2Ct . Reactor antineutrino limit is from 
Ref. (Hj. The limit obtained by CHS H2 is also shown. 



One of the open questions in neutrino physics is under- 
standing the role of mixing between the first and third 
flavor generations, #13. In this regard we also explored 
if the uncertainties coming from the lack of knowledge 
of and the counter-term Lia are comparable. In the 
limiting case of small cos #13 and Sm^i 3> Sm^i , which 
seems to be satisfied by the measured neutrino proper- 
ties, it is possible to incorporate the effects of 6*13 rather 
easily. In this limit the three-flavor survival probability 
is is given by H 0111 



f3x3(fe ->■ V & ) = COS* $13 Plyi.l{v e ~> V e with N e COS 2 9 13) 

+ sin 4 6»i 3 . (8) 

In Eq. (|SJ the quantity P2x2(i / e v e with N e cos 2 #13) 
is the standard two-flavor survival probability calculated 
with the modified electron density N e cos 2 6*13 and the 
standard initial conditions. This suggests that for small 
values of $13 the survival probability and consequently 
the counting rate can be linearized in cos 4 #13: 



Count Rate ~ A + B (1 - cos 4 9 13). 



(9) 



The neutral- and charged-current counting rates lin- 
early depend on L±a while elastic scattering rate does 



not. Conversely the charged-current and elastic scatter- 
ing rates linearly depend on cos 4 #13 while the neutral- 
current rate does not. Hence it is reasonable to compare 
their relative contributions. To this end in Figure r4] we 
show the allowed 6*13 and Lia parameter space when #12 
and 5m\ 2 are taken to give the minimum % 2 values to 
reproduce the data. Results where the fixed values of #12 
and 8m\ 2 obtained using only the best fit of the SNO 
data (left-hand panel) and using the best fit of all the 
solar neutrino data along with KamLAND results (right- 
hand panel) are both shown. The dark-shaded region 
corresponds to the case when Jb — 1. The light-shaded 
region corresponds to the case when the 8 B flux is uncon- 
strained. We observe that the uncertainty coming from 
the lack of a better knowledge of #13 is larger than un- 
certainty coming from not knowing Lia precisely. One 
also observes that as 6*13 increases the allowed L\a re- 
gion shifts toward larger values of L\a- When 9 13 = 
the electron neutrino flux is lost into only one channel: a 
particular linear combination of /i and t neutrinos |36j| . 
But when #13 7^ additional flux is lost into the or- 
thogonal channel as well. The slightly decreased electron 
neutrino survival probability reduces the charged-current 
and the elastic scattering count rates so that a slightly 
larger L\a is needed to compensate the resulting decrease 
in the count rates at each bin. 




FIG. 4: Allowed parameter space when 9i2 and 5rn\ 2 are fixed to give the minimum \ 2 values to reproduce the SNO day-night 
spectrum (left-hand panel) and all solar neutrino experiments along with the KamLAND experiment (right-hand panel). The 
shaded areas are the 90 % confidence level region. 95 % (solid line), 99 % (log-dashed line), and 99.73 % (dotted-line) confidence 
levels are also shown. The dark-shaded region corresponds to the case when the 8 B flux is fixed to be the Standard Solar Model 
value. The light-shaded region corresponds to the case when that flux is unconstrained. 



In conclusion we showed that the SNO experiment with 
increased statistics using additional input from other so- 
lar neutrino experiments can significantly reduce the un- 
certainty in determining the precise form of the axial two- 
body currents at low energies. We also showed the con- 
tribution of the uncertainty in L\a to the analysis and 
interpretation of the SNO data is nearly negligible. The 
effect of this uncertainty is smaller than effects of a non- 
zero value of 6*i3 or even than effects of possible solar 
density fluctuations |33|. Finally our most conservative 
value for L\a is significantly larger than that was ob- 
tained by CHR. One reason for this may be the treatment 
of neutral- and charged-current count rates together in 
the global analysis. 
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